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A narrow state, which we label D sJ (2458) + , with mass 2458.0 ± 1.0 (stat.) ± 1.0 (syst.) MeV/c 2 , 
is observed in the inclusive Dj~7r°7 mass distribution in 91 fb _1 of e + e~ annihilation data recorded 
by the BABAR detector at the PEP-II asymmetric-energy e + e~ storage ring. The observed width is 
consistent with the experimental resolution. The data favor decay through Dg (2112) + 7r rather than 
through D*j(2317) + 7. An analysis of D+ir° data accounting for the influence of the D 3 j(2458) + 
produces a D*j(2317) + mass of 2317.3 ± 0.4 (stat.) ± 0.8 (syst.) MeV/c 2 . 

PACS numbers: 14.40.Lb, 13.25. Ft, 12.40.Yx 
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Interest in the spectrum of charmed mesons has been 
heightened by the discovery by this collaboration p| of 
a narrow state, produced in e + e~ — > cc collisions at 
the PEP-II collider, decaying to D+n" 2], with mass 
2317 MeV/c 2 , approximately 41 MeV/c 2 below the DK 
mass threshold. This state, the -D*t(2317) + , has been 
confirmed by CLEO Q and Belle 0, ^. Along with 
the D*j(2317) + , we noted [| the presence of a narrow 
peak in the D+ir°j mass distribution near 2.46 GeV/c 2 . 
Because this signal is near the kinematic overlap of the 
L>* J (2317)+7 and D*(2112)+7r° systems, special atten- 
tion is required to remove associated backgrounds and to 
distinguish between the two possible decay modes. Such 
an analysis is the subject of this paper. 

This state near 2.46 GeV/c 2 has been seen by CLEO Q 
and Belle Pjj in the inclusive Z?+7r°7 mass spectrum and 
by Belle pj in exclusively reconstructed B decays. 

To investigate the D+7r°7 spectrum, we study Df can- 
didates from e + e~ — > cc (at a center-of-mass energy near 
10.6 GeV) that decay to K~ K + tt + . Particle identifica- 
tion is used to provide clean samples of charged K and 
7r candidates, which are combined using a geometric fit 
to a common vertex. Backgrounds are suppressed by se- 
lecting decays to K*°K + and 4>ir + . A description of this 
sample and additional details can be found elsewhere . 
Events with 1.954 < m(K-K + Tr+) < 1.981 GcV/c 2 are 
taken as D+ candidates. 

A candidate 7r° is formed by constraining a pair of 
photons each with energy greater than 100 MeV to em- 
anate from the intersection of the Df trajectory with the 
beam envelope, performing a one-constraint fit to the ir° 
mass, and requiring a fit probability greater than 5%. A 
given event may yield several acceptable 7r° candidates. 
We retain only those candidates for which neither photon 
belongs to another otherwise acceptable 7r°. 

Each candidate is combined with all combinations 
of accompanying ir° candidates with momentum greater 
than 300 MeV/ c and photon candidates of energy greater 
than 100 MeV. To suppress background, photons that be- 
long to any ir° candidate are excluded and we require the 
momentum, p*, of each Dfir 0r y combination in the e + e _ 
center-of-mass frame to be greater than 3.5 GeV/c. The 
last requirement also removes any D+7r°7 combination 
from B decay. 



The Djn 7 invariant mass distribution is shown 
in Fig. Ui. A clear enhancement is observed near 
2.46 GeV/c 2 . The background underneath this peak is 
from several sources, which can be described in terms of 
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FIG. 1: (a) The mass distribution for all selected D£ 
combinations. The shaded region is from Df sidebands de- 
fined by 1.912 < m{K~K + -K + ) < 1.933 GeV/c 2 and 1.999 < 
m(K-K + Ti + ) < 2.020 GeV/c 2 . (b) The value of Am 7 ver- 
sus Am n o for all combinations. The horizontal lines delineate 
three ranges in Am 7 . (c) The Am,o mass distribution for the 
middle range of Am 7 (points) and for the average of the up- 
per and lower ranges (shaded histogram), (d) The difference 
between the two distributions shown in (c). The curve is the 
fit described in the text. 



mass differences defined as 



Ato 7 = m{Dg 7) — m(Dg) 



Am-o 



mp+ 7 7r°) - m{D+i) 



(1) 
(2) 



A scatter plot of the data is shown in Fig. ^3. Particular 
background patterns are visible: D*(2112) + — » Df^y de- 
cay combined with an unassociated ir° , which appears as 
a horizontal band, and Z?*j(2317) + — > Dfir decay com- 
bined with an unassociated 7, which appears as a band 
that is almost vertical. 

To demonstrate the existence of a signal above these 
backgrounds, the upper histogram of Fig. shows 
D+7r°7 combinations in the D*(2112) + signal region, and 
the gray histogram, scaled to the area of the signal re- 
gion, corresponds to the two Z?*(2112) + sidebands. We 
conclude that a signal for a state decaying to Dfir°j ex- 
ists over a background resulting from D*j(2317) + and 
an unassociated 7. This background peaks at a mass 
slightly higher than that of the signal. A Gaussian fit 
to the subtracted mass distribution (Fig. |TJI) indicates a 
narrow signal at Ato^o = 346.2 ±0.9 MeV/c 2 (statistical 
error only). 

The state corresponding to this signal, which we label 
L> sJ (2458)+, may decay to £>+tt°7 through D*(2112)+7r° 
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or D*j(2317) + 7. To disentangle these modes and reli- 
ably extract the parameters of the signal, we apply an 
unbinned maximum likelihood fit simultaneously to the 
£)+7r°7, -D+7T , and Df-f invariant masses of all Z?+7r°7 
combinations using the channel likelihood method 0. 
This fit describes the probability density function of the 
two D s j(2458) + decay channels as the product of a Gaus- 
sian shape in the D+ir°j mass distribution and a Gaus- 
sian shape projected into the Dfn a or 7 mass axes, 
as appropriate. Because the daughter resonances are 
narrow, interference between the two D s j (2458) + decay 
modes cannot be resolved, and so is ignored. 

Sources of background in the -D+7r°7 spectrum in- 
cluded in the fit are purely combinatorial background 
(Df meson combined with an unassociated ir° and 7), 
£>*(2112) + — ■> Dfj decay combined with an unasso- 
ciated 7T°, and L>* 7 (2317) + — > Dfn° decay combined 
with an unassociated 7. The fit also includes a con- 
tribution from D s ,/(2458)+ -> D*(2112)+7r° decay but 
with an unassociated 7 replacing the 7 from D*(2112) + 
decay. The fit determines the relative size of the back- 
ground and signal contributions, the mass and width of 
the £> sJ (2458)+, and the D* 7 (2317)+ mass. 

The likelihood fit is validated using Monte Carlo (MC) 
simulation. This simulation includes e + e~ — > cc events 
and all known charm states and decays, including the 
D*j(2317) + and the signal under study. The gener- 
ated events were processed by a detailed detector sim- 
ulation and subjected to the same reconstruction and 
event-selection procedure as the data. 

As shown in Fig. the fit provides a good descrip- 
tion of the Dfir°j mass distribution observed in the 
data. The D s j (2458) + signal for a particular decay mode 
can be isolated by calculating a weight for each D+n°~f 
combination proportional to the relative likelihood con- 
tributed by the decay mode of interest. Distributions 
of events so weighted can be compared to the likelihood 
function to validate the fit. This is shown in Figs. |2b 
and 13;. A \ 2 probability calculation gives 22%, 74% and 
11% for fig. [2]a, b and c respectively. The resulting yield 
of correctly reconstructed Z? sJ (2458) + -> £>*(2112)+7r° 
(£> sJ (2458)+ -► D* J (2317)+7) decays is 195±26 (0±23), 
consistent with the fit shown in Fig. [IJi. Excluding the 
D s j(2458) + from the likelihood fit decreases the log like- 
lihood by approximately 57, corresponding to a signifi- 
cance of more than 10 standard deviations. The fit yields 
a D 5 ,/(2458)+ mass of 2458.0 ± 1.0 MeV/c 2 with an rms 
width of 8.5 ±1.0 MeV/c 2 . 

The likelihood fit uses the shapes of the D+ir° and 
D^"7 mass distributions to distinguish between the two 
possible decay modes, L» sJ (2458)+ -> D;(2112)+7r° and 
LW(2458)+ -> L»; j (2317)+7. These shapes are influ- 
enced by the kinematic constraints of D s j(2458) + de- 
cay shown in Fig. Figs. [3jD~|2t; show the sideband- 
subtracted Dfn and -D+7 mass projections compared 
with MC simulations of the two hypotheses (scaled to 
match the data yield). The L> sJ (2458)+ -> L>;(2112)+7r° 
decay mode (solid histograms) produces a narrow Df"f 
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FIG. 2: Maximum likelihood fit results overlaid on the 
7?j"7r°7 mass distribution with (a) no weights and after apply- 
ing weights corresponding to (b) the decay D* (2112) + 7r° and 
(c) the decay D*j (2317) + 7. (d) The mass spectrum of D^n° 
combinations (with no 7 requirement). The solid curve is the 
fit described in the text. The dashed and lower solid curves 
are the contributions from D s j (2458) + decays and combina- 
torial background, respectively. 



mass distribution and a wide Dfir mass distribution. In 
contrast, the D s . / (2458)+ -> D* J (2317)+7 decay mode 
(dashed histograms) produces a wide Dfj mass distri- 
bution and a narrow Dfir mass distribution. Figures^ 
andEt show that the £> sJ (2458) + -> D*(2112)+7r° hy- 
pothesis is in better agreement with the data. 

Our previous measurement |]J of the D*j(2317) + mass 
using the decay £>*j(2317)+ -> D+tt° did not explicitly 
consider background from D s ,/(2458)+ -> £>*(2112)+7r° 
decay. This background peaks in the Dfir mass spec- 
trum just below the D*j(2317)+ mass. Shown in Fig. 01 
is the Dfn invariant mass distribution for a sam- 
ple of -D+ candidates combined with all ir° candidates, 
with p* > 3.5 GeV/c. Superimposed on this distri- 
bution is a binned fit that includes the contribution 
from the D s j(2458) + as estimated from MC simulation 
and a quadratic background function. The result is a 
D* J (2317)+ yield of 1022 ± 50 events, a mass of 2317.3 ± 
0.4 MeV/c 2 , and measured rms width 7.3 ± 0.2 MeV/c 2 . 
These results are an improvement over our earlier mea- 
surement 0. 

We divide the sources of systematic uncertainty in the 
Z? s j(2458) + and £>*j(2317) + mass values and production 
rates into three categories. The first category is associ- 
ated with the fit procedure. Likelihood fits to MC sam- 
ples that include samples of £> sJ (2458)+ -> £>*(2112) + 7r° 
and £> s j(2458) + — > £>*j(2317) + 7 decays correctly repro- 
duce the given sample sizes within statistical errors. The 
average values of the fit results obtained using statisti- 
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FIG. 3: (a) The Df^y versus Dfir mass distribu- 
tion for all D^~7r°7 combinations. The decay of a zero- 
width D s j(2458) + is kinematically restricted to the re- 
gion between the two curves. (b) Sideband-subtracted 
D])~7 mass distribution with MC simulation for (solid his- 
togram) D sJ (2458) + -> D*(2112) + 7r° and (dashed his- 
togram) D sJ (2458) + -> D*j(2317) + 7. (c) A similar plot for 
the Dfir mass distribution. 



cally distinct MC samples corresponding to the measure- 
ments in the data are used to place limits on any fit bias. 

We obtain the background distribution used in the like- 
lihood from a random selection of Df , 7r°, and 7 candi- 
dates taken from the MC D+tt°j sample. To test our 
sensitivity to this distribution, various selection require- 
ments are altered within reasonable bounds to provide 
alternate background samples for use in the fit. The re- 
sulting changes in yield and mass are used as the second 
category of systematic uncertainty. 

Reconstruction of the decay sequences is the third 
source of systematic uncertainty. To evaluate the reliabil- 
ity of the MC determination of 7r° efficiency and momen- 
tum calibration, we use control samples of K$ — ^ 7r°7r° 
and r — > tt°X. On this basis, we assign a systematic 
uncertainty of ±5% in ir° reconstruction efficiency and a 
relative ±1% in 7r° momentum bias. Similar studies for 
7 reconstruction reveal a systematic uncertainty of ±3% 
in 7 reconstruction efficiency and ±1% in energy bias. 
Uncertainties in the Dj~ and £)*(2112) + masses, taken 
from world averages fgj], also contribute to the system- 
atic uncertainty. 

The resulting total systematic uncertainty in the 
L> sJ (2458)+ [Z»* i/ (2317)+] mass is ±1.0 [±0.8] MeV/c 2 . 

Using the yields from our fit and correcting for effi- 
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ciency, we estimate the relative production rate 

= a(^, 7 (2458)+)g( J P sJ (2458)+ -> j;(2112)+7r°) 
( 7(£>: J (2317)+)B( J D: J (2317)+ - D+ifi) 

(3) 

to be 0.25 ± 0.03 (stat.) ± 0.03 (syst.), requiring p* > 
3.5 GeV/c for both states. We also estimate, at 95% 
C.L, 

g(J> aJ (2458)+ D; j (2317)+ 7 ) 
B(Z>.j(2458)+ -f D;(2112)+7T°) ' 1 ' 



The observed rms width of the Z? s j (2458) + is consis- 
tent with detector resolution, as determined by Monte 
Carlo studies. We conclude that the intrinsic width of 
the L> sJ (2458)+ is small (r < 10 MeV/c 2 ). 

The mass of the D s j (2458) + lies above DK and be- 
low D*K thresholds. The narrow width and the isospin- 
violating decay to Z?*(2112) + 7r° indicate that decay to 
DK is forbidden and suggest an unnatural spin-parity 
assignment for the state. Belle has observed the decay 
-D s ,/(2458) + — > -D+7 introduction from both cc contin- 
uum £| and B decay [5J. Such a decay rules out J = 
and favors a 1 + interpretation. Decay distributions stud- 
ied by Belle further support J = 1 for D s j (2458) + and 
also J p = 0+ for D* J (2317)+. The apparent absence 
of the decay L> s j(2458)+ -> D*j(2317) + ^ may indicate 
that the electromagnetic decay mechanism cannot com- 
pete with D*(2112) + 7T°, which may be a strong, but 
isospin-violating, process resulting from r]-ir mixing, as 
discussed by Cho and Wise [j| . 

Our measurement of the -D s j(2458) + mass (2458.0 ± 
1.4 MeV/c 2 , with combined statistical and systematic un- 
certainties) agrees with that obtained by Belle (2456.5 ± 
1.7 MeV/c 2 ) 3, but is two standard deviations smaller 
than that obtained by CLEO (2463.1 ± 2.1 MeV/c 2 ) 0. 
We obtain a relative yield (R — 0.25 ± 0.04) which agrees 
with that of Belle (0.26 ± 0.08). Both values are some- 
what smaller than that reported by CLEO (0.44 ±0.13). 
Our reanalysis of the D*j(2317) + — ► Dfir sample to ac- 
count for background from the D s j(2458) + gives a mass 
of 2317.3±0.4 (stat.)±0.8 (syst.) MeV/c 2 , which remains 
consistent with results from CLEO Q and Belle p|. 

In summary, in 91 fb _1 of data collected from the 
BABAR experiment, we have observed a narrow state that 
decays to Dfir ^ with a mass of 2458.0 ±1.0 (stat.) ± 
1.0 (syst.) MeV/c 2 . The only significant D+7r°7 decay 
mode we observe is through D*(2112) + 7r°. We measure 
a mass and yield relative to the Z)* 7 (2317) + similar to 
those measured by Belle though smaller than those re- 
ported by CLEO. The observed width is compatible with 
our mass resolution. After including the influence of this 
state, our new measurement of the -D*j(2317) + mass is 
2317.3 ± 0.4 (stat.) ± 0.8 (syst.) MeV/c 2 . 
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